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SURMAEY 


A strairt-energy theory is developed for the calculation of the 
buckling load in pure bending of reinforced monococLue cylinders wMch 
have a symmetric cutout on the conpression side and buckle according to 
the general instability pattern. Coii 5 )ufcaticKis are carried cut for the 
cylinders tested earlier at the Polytechnic Institute of Brooklyn 
Aeronautical Laboratories . The theoretical curve is similar in shape to 
that obtained experimentally, but the theoretical values are consistently 
too high. The .deviation ia 39-3 percent in the worst case. 


IHTEODUCTION 

General" instability is defined as "the simultaneous buckling .of . the. 
lon^tudinal and Ol'rcumf erant ial reinforoing elements .of armonocoque. 
cylinder together with the sheet attSiChed to th^. The general instabil- 
ity of reinforced circular monocoque cylinders, subjected to pure bending 
has been investigated in some detail at Polytechnic Institute of Brooklyn 
Aeronautical Laboratories’ and Guggenheim. A.erpnautlcal Laboratory, 
California Institute of Technology, under the sponsorship of the National 
Advisory Committee for Aeronautics (references 1 to -8) . This theoretical 
and exp er imental work dealt with coaplete cylinders not ha'.’-ing cutouts - ’ 
It can be expected. that a cutout decreases the buckling load in general 
instability since part of the elastic support is lost , when a portion of 
the structure is removed. This conjectxure was verified in recent exper- 
iments carried out at Polytechnic Institute of “Brooklyn Aeronautical 
Laboratories which dealt with the general instability of and the stress 
distribution in monocoque cylinders with a syimaetri.o cutout. Eeference 
9 contains a report on these experiments. A theoretical study of the- 
stress distribution in the cylinders is presented in reference 10. 
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Beference 11 deals with an experimentetl investigation of cylinders having 
a side cutout. 

In the present report the huckling load of reinforced monocoque cyl- 
inders with a symmetric cutout on the congpression side is calculated by 
strain— energy methods. The deflected shape at buckling is represented 
by a full sine wave extending over the length of the cutout in the axial 
direction and by the first seven teims of a Fourier expansion in the cir- 
cumferential direction. The circumferential coordinate is measured from 
the edge of the cutout and the length of the interval in which the 
Fourier series is defined is considered as one of the parameters of the 
problem. The boundary conditions at the end of the Interval determine 
four of the seven coefficients of the series while one of them is inde- 
terminate as in all bucklii:g problems. The remaining two coefficients^ 
as well as the wave-length parameter, are calculated frcm the requirement 
that the buckling load be a minimum. 

The following strain-energy quantities are considered: radial and 

tangential bending as well as torsion of the stringers; bending of the 
rings in their plane; and shear in the sheet. The extensional strain 
energy stored in the sheet is taken into account by adding an effective 
width of sheet to the stringers and the rings. In the calculation of 
the work of the external forces a linear force distribution is asstuned 
in preference to a lineair strain distribution in bending. This asexmop— 
tion is in better agreement with the experiments described In reference 9 * 

The buckling load is calculated from the requirement that the strain 
energy corresponding to the transition frcm the unbuckled into the buck— 
le'd shape be equal to the work done by the applied loads. The minimum 
vaJLue of the buckling load is found by assuming the circumferential wave 
length to be equal to the length of eome integral number of stringer 
fields, calculating the values of the two Fourier coefficients that min- 
imize the buckling load in the case of the assumed wave length, determin- 
ing the buckling load, and conparlng it with values obtained on the as— 
Buugptlon of other different wave lengths. The final results of the 
numerical work are presented in the fom of buckling loads calculated for 
three different circumferential wave lengths for each of the three sizes 
of the cutouts tested. In each case these buckling loads define a mini- 
mum. All the calculations were carried out for one— half the cylinder 
because of the symmetry of both structure and loading. 

For a Substantial share in the numerical work the authors are in- 
debted to Bernard Levine. The Investigation was conducted under the 
spotieorshlp and with the financial aid of the National Advisory Committee 
for Aeronautics. 
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SYMBOLS 


a,ao,ai,a2,aa 

A 


'bjlJij'baj'b-a 

C 

d 


E 


G 


^eff 

i 


I 


Ir 


^strj. 


^strt 

J 


ki,lC2,]s3,l£: 

L 

Li 

m 

M . 


Fourier coefflclexits 

cross-sectional area of a stringer plus its effective vidtli 
Fourier coefficients 

geometric factor in toiBional rigidity GC 
width, of panel measured along the circumference 
Young's modulus 
shear modulus 

shear modulus of sheet covering at zero compressive load 
effective shear modulus 

index indicating position along circiunference 
moment of inertia 

moment of inertia of ring section and its effective width 
of sheet for tending in its own plane 

moment of inertia of stringer section and its effective 
width of sheet for tending in the radial direction 
(about a tangentleil axis) 

moment of inertia of stringer section and its effective 
width of sheet for tending in the tangential direction 
(about a radial axis) 

index indicating position along axial direction 

trigonometric functions of cp, x, n, a, t 

length of cutout 

distance between adjacent rings 

number of rings involved in the failure 

applied bending moment; function of n, a, b appearing in 
the strain, energy of bending in the rings 
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Mcr applied TDending moment at general instability 

n parameter defining the length, of the general instability 

in the circumferential direction 


N = 0 .0275 [(2itr/d) + l] 


Pi^ P2 
^cr 


Q 


r 

E 


s 


S 

t 

U 

Ur 


Uah 


^strt 

Ut 

2^f 


polyaomlal functions of a and b 

force ceirried by one of the stringers at the edge of the 
cutout at general instability 

force carried by the ith stringer 

function of x appearing in the shear strain energy 

radius of cylinder 

function of • 9 , n, a, b appearing in the shear stn.in 
energy 

number of stringers involved in one -half Uie general- 
instability bulge 

total number of stringers in the cylinder 
thicJmeae of sheet covering 
strain energy 

strain energy stored in tho rings because of bending (of the 
rings) in their o'wn plane 

strain energy stored, in the sheet covering because of 
shear 

strajji energy stored in tho stringers because of bending 
about a tangential axis 

strain enei'gy stored in the stringers because of bending 
about a radial axis 

strain energy stored in the stringers because of torsion 
effective width of sheet 

radial displacement of a point on a ring or a stringer 
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tangential olsplacemsiit- of a point on a ring or a stringer 
•fforic dosae by the applied forces 

coordinate measuring distance along the axis of the cylin- 
der from the edge of the cutout 

cutout angle 

coefficients used, in the calculation of the shear strain 
in a panel due to displacements of its comers 

shear strain 

distance of neutral axis from horizontal diameter of cyl- 
inder 

normal strain in a stringer 

huclcling strain of a panel of sheet covering 

angular coordinate vith origin at the edge of the cutout 

tee lEBlSCTED SHAPE 


In the eaperiments described in reference 9 it was observed that at 
buckling the ware length in the axial direction was almost exactly eq.ua! 
to the length of the cutout . For this reason it is assumed In the theoay 
that the rings .bordering the cutout are rigid in their planes. The cylin- 
der is then thought of as being cut through these rings end the external 
moments are applied in the sections. ¥ith the notation of figure 1 the 
distorted shape of the strirgere le assumed to be 

Wr = ao(kx/2)[l — oca (2tcx/l) ] = aoki sinS (jtx/L) Cl) 

where Wj* is the radial deflection, and ki a prcrportionallty factor ■ 
dependent upon the angle cp.. 

The olrctonferentlal wave' pattern could not be determined trith suffi- 
cient accuracy in the tests. It is assumed, therefore, to be represented 
by the foUcwing trigoncanetric erpression: 

Wr = ka[ (ao + a^ cos n£p + ag cos + as cos + bi sin nep 


+ hs Bin 2rtcp + ha sin 3n»p) I (2) 
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wliere . ■ JLs a proportionality faojbor .dependent upon . x. (Because of 

eq.uatlon (l) kg is ein.2 (jti/L).) Equation (2) is valid, provided 

0 < cp < rt/n * (2a) 

When cp is greater than ■ n/n, the deflections are assumed to he zero. 
Consequently n is the parameter defining the length of the bulge. 

Since in thin rings extensional deformations Involve much more 
strain energy than do bending deformations, the deflections of the rings 
are assumed to be inextensional. This assurgptlcn determines the tangen- 
tial displacements wt when the radial displacements are given. The 
connection between the two was developed in reference 3 aJid stated in 
equation ( 4 a) of that reference: 

Wr = (3) 

It follows from equations ( 2 ) and ( 3 ) that w-^ may be taken as 
wt = kaC-sk)? ~ (a-i/n) sin nrp — (aa/2n) sin 2ncp — (ag/3n) sin 3 i^ 

+ (bi/n) cos ncp + (b2/2n) cos 2ncp + (bs/3u) cos 3 ucp ] ( 4 ) 


provided 


0 < qp < jt/n 

Because of the symmetry cf both structure and loading these expres- 
sions are equally applicable when the angle (p is measured from either 
one, of the edge stringers. An obvious limitation of the formulae is 

V i 

■ . . , a + < It. (5) 

If it is required that there be a smooth transition between the 
bulge and the nondistorted port of the cylinder at tp =» .{s'/n), then the 
followiig conditions must be satisfied: 

(1) The tangential displacement must vanish: 

. wt = p . . when cp *=.jt/n . . . , . ^ : ( 6 a) 

(2) The radial displacement imist vanish: • 

Wj. =.0 . when cp =? it/n . .. . - ( 6 b) 

(3) There must be no. sudden change in the direction of the tangent: 
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Swr/09 =» 0 wiieii tp = ff/n • 

(W) There must he no sudden change in the curvature: 
S Wp/Scp^ = 0 when 9 = rt/n 


(6c) 

(6d) 


The mathematlcel foraulation of these requirements was discussed in 
detail on pages 10 and 11 of reference 3- four conditions contained 

in aquations (6a) to (6d) estahlish four relationships between the 
Fourier coefficients in equation (h) and make it possible to express any 
four coefficients by means of the remaining three. If ao, and bi 
are retained as the basic parameters, the following foiir equations are 
obtained; 

&2 = (8/5)c.i - ( 9 / 5 )ao 
as = (3/5) ax — (V 5 )ao 

> ( 7 ) 

ba * (l6/p)bx + (18/5 )nac 
ba = ( 9/5 )tx + (l2/5)nao 

With the notation 

(ax/ao) = a and (bx/ao) = b (8) 

and after substitution of the expressions contained in equations (7), 
a combination of equations (l) and (2) gives for the radial displacement 

wr = aokx sin2 (nx/L) (9) 

where 


hx = [1 + a cos ncp + (1.6a — 1.8) cos 2n9' 

+ (0.6a — 0.8) cos 3p9 + b sin xxp + (3.2b'+ 3-6rt) sin P-ntp 
+ (1.8b + 2. hit) sin 3 tKp ] (S®-) 

Similarly the tangential displacement becomes 

w^ = sin^ (ioc/L) (lO) 

ks = (l/n) [-09 — a sin n9 — (I/2) (l.6a — 1,8) sin 2n9 

~ (l/3)(0*6a — 0.8) sin 3^9 + b cos 339 + (l/2)(3.2b + 3-6it) cos 2n9 
+ ( 1 /^( 1 . 8b + 2. hit) cos 3n9 3 (lOa) 
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Eq,uations (9) and (10) Are valid, provided 

0 < cp < n/n (10b) 

When qp is greater than n/n the deflections are asstaned to vanish. A 
typical example of the deflections at buckling in the plane of the rings 
is shown in figure 2. 


CALCULATION OF TBE STEAIN ENEEGY 
Strain Energy Stored in the Eings 


The strain energy stored in half of any ring is 


U = (1/2) (El)r/r 


(El)r/r® 


Vr + (S^Wj^/dcp^) 


dxp 


( 11 ) 


in accordance with eq.uations (c) on page 11 and (7) page 12 of refer- 
ence 3* Substitution of the value of Wr from. eg.uation (9) and summa- 
tion over all the rings contained in the axial wave length yield 


Ur = (1/2) (ao^/r®) ) (El)r sln^ 


(^Jtxj/lj ^ [1 + a(l — n®) cos ncp 
+ (1.6a — 1.8) (1 — kn^) cos 2ncp + (0.6a. — 0.8) (1 — 9a^) cos 5ikp 


+ b(l — n^) sin n<p + (3-2b + 3»6jt)(1 — 4n^) sin 2ncp 
+ (1.8b + 2.4jt)(1 - 9n^) sin 3iifp3 ^ (12) 


where Ur is the strain energy stored in all the rings in one-half of 
the cylinder. The subscript j refers to the individual rings the 
total number of which is m within the lei^gth of the cutout. If the 
integration is carried out and the value of the definite integral is de- 
noted by M it is possible to write 
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nM = [jt + 10.053096(1 - 9nS) + 206.01005(1 - 4n2)2 + 90.30338 t (1 ~ 

- 18.095573 (1 - (1 - 9n^) ] 

+ a [-9.0477868(1 - 4n2)® - 1.507^45(1 - 9n-)^ 

+ 30.159289(1 - a^) (1 - 4a^) -i 18.095573 (1 - 4a^) (l - 9n^) ] 

+ b [4(1 - n ^) + 2.4 0. - 9n^) + 113.69784 (1 - 4a^)^ + 42.636690(1 - 9n^) 

+ 2.4(1 - n®) (1 - 4n®) - 3.68(1 - 4n'^) (1 - 9n^)-] 

+ a2 [(V2)(l - n2)2 -+ 4.0212386(1 - 4n2)a + 0.5654867(1 - 9D .^ f ] 

+ [(jt/2) (1 - a^) ® + 16.084954(1 - 4a®) ® + 5.08938(1 - 9a®)®] 

+ ab [6.4(1 — a®) (1 — 4n^) + 3-84(1 — 4n^) (1 — 9a^) ] (13) 


The strain energy Is therefore 

m. 

Uy = (1/2) (ao®/r'^)M / (El)^ sln^ («3 :,/l) (l4) 

J=^ 

When the bending rigidity (El)x* is the same fer all the rings, the 
summation yields a result in closed form as vas shown in the appendix 
of reference 12 : 

m 

^ sin"* («xj/L) = ( 3 / 8 ) (m + l) (15) 

J=i 


provided 


m > 1 


(15a) 


t/hen m = 1 the value of the Bxuinnatlon is 1. The strain energy of 
bending stored in all the rings is consequently 


Ur = (3/16) (ao^r"^) (El)r(m + l)M (16) 


In eqmtion (I 6 ) the value of M d^en^ upon n, a, and b. Values of 
M ccn^iuted for n = 4, 2 . 666 . . ., 2 , 1 . 6 , 1.333. - . cKirrespcsiding to 
s = 2, 3j 4, 5 , 6 , respectively, are listed in table I. These values of 
n correspond to buckling patterns in which the bulge ends at one of the 
16 Btringera contained in the spectmsns tested. 
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Strain Energy Stored in the Stringers 

The strain energy stored in the stringers because of bending in the 
radial direction is 

L 

^ 7 > s> 2 

Ustrj. = ^ (1/2) (El)atrj. J w^/Sx ) dx (IT) 

where the summation is extended over all the stringers contained in one~ 
half the cylinder. Substitution of the value of Wj. from eq.uation ( 9 ) 
into eq^uation (IT) and integration yield 

L 

Ustrr = ) (l/2)(El)atrr f ^ fci®(2it^L®) ^ cos^ (2nxA) 

£-J Jq 

= ao^ (n /L ) V' (El)etrr (^) 


The moment of inertia letrr s'trlnger varies around the 

circumference of the cylinder because the effective width of the sheet 
to be added to the stringer section changes . The vaduee of Istrr vere 

determined, for each of the cutout sizes investigated, according to the 
principles stated in reference 1. Similarly ki^ was confuted for each 
stronger. 


The strain energy stored in the stringers because of bending in the 
tangential direction is 


L 

Ustrt = ) (l/2)(EI)strt / (&®wt/Sx^) dx 

' 'o 


( 19 ) 


where the summation is extended over all the stringers contained in one- 
half the cylinder. With the aid of equation (10) the sbrain energy can 
be given as 


Ustrt = ao®(n‘‘/L'') Y t3^(El)strt (20) 

I I 

Since both kg and Istrt vary from stringer to stringer the summation 

indicated in equation (20) was evaluated numerically for each cutout size 
investigated. 


The strain energy stored in the stringer because of torsion is 
^ L 

^t (l/2)GC ^ [(SSrr)/(^x dcp) ] dx 


( 21 ) 
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In this equation' (l/r) 5 <p) i's tiie unit ’an'gi© of twiat of 

the stringer, and the summation. imist he carried oht- oyer ell the 
"stringers contained in one-half "the cylinder. Iii "the eacpression for 
"the Saint -Venant torsional rigidlty,- 

C = 0 .l 4 a* ( 21 a) 

since the test specimens were provided "Wi"th. sqtLare section stringers 
of edge length a. Differentiation gives 

(S%y)/(3x Sep) = ao3s:4(jt/L) sin (2rtx/L) ^ (22) 

where . " ' 

IC4 = n[j- a sin mp - (3 .2a - 3.6) sin 2n<p - (1.8a - 2 ,k) sin 3^^ 

+ h cos nfp 4 ( 6 . 4 h +■ ?. 2 ic) cos 2 nq> 4 ( 5 »‘**^h + T* 2 jc) cos 3^19]' ( 22 a) 


Hence the strain, energy of torsion is - - ■ • 

Ut = [(GC)/(Lr 2 )J . (^) 

where the summation is extended to include all the stringers contained 
in one-half "the cylinder, Eie term GC is hefore the summation sign 
in equation (23) since according "to Sai.nt "Venant the variation of the 
torsional rigidity, caused hy the different amounts of effective 
■Hld"fch of sheet, is so small that it "was considered permisslhle "to 
assume GC a cons"tant. Again a numerical evaluation of the stunraa- 
tion was carried out "when the' s"train energy was calculated . 


Strain Energy of Shear Stored in the Sheet 

The shear s"fcrain energy in a panel is taken as heing proportional 
"to "the average effective shear modulUB <j0f f multiplied hy the square 

of the average shear strain 7 in the panel, Ihe latter is-. calcu- 
lated from the dlsplacemen"bs of the four comers of the panel as was 
done in reference 3 . Then "the total strain energy of shear stored 
in "the sheet is . • • 


Ush = (172)^^72 Goff td ( 2 k) 

vhere "the summation extends over all the panels con"tained in one- 
half "the cylinder. ■ - 
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The value of the effective shear modulus depends upon the matjnitude 
of the con 5 >r©ssive strain in tho panel, as was shown in reference 13- 
The eiiplrlcal formula recommended there for the computatlcn of Gepf is 

Geff/O-o = (1 ~ + N (2p) 

where 

N = 0.02T5E2n(r/d) +1]' (25a) 


and Go is the shear modulus in the absence of compressive stresses, 
e the con^jressive strain prevailing in the panel, and the com- 
pressive strain when the panel of sheet Duc^ee. 


Since the displacements in the axial direction are small and of the 
second order, the displacements of the comers of the panel need he in- 
vestigated only in the plane of the rings. The effect of rotation of 
the ring upon the shear was neglected. Portaulae for the calculation of 
the shear strain from the displace^nts of the comers were developed in 
reference lO and were' presented in figure 23 of that reference. With 
the notation and sign convention of figure 3 of the present report the 
shear strain is 


r - (orAi) jvri^ J - J 


where the first suhscript refers to the circumferential, and the second 
to the axial location of the comer of the panel. In reference l4 the 
values of the numerical factors Oif and at were determined. Simpli- 
fied formulas were given on page 27 of reference l4 which represent 
these factors very accurately when the angle ‘ d/r is of the order of 
magnitude found in monocoque fuselages. The formulas can 'be written in 
the following slightly changed form: 


Or = ( 1 / 10 ) (d/r) = (1/10) (2«/s) 

CH = ~ 1/2 


(27) 




L i 

j=o i=o 


(^eff/^o^i 


Substitutions yield 

Ugi, = (l/2)(td/Li)Go 


(28) 
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where Q is a function of i only, and E a function of cp only. The 

suffluaticn ^ Q gives a result in closed form: 
m ^ 

V Q* = y ‘fsih^ [(«j)/(m + 1)] — sin® [x (j + l)/m + l) l^i- 

= (lA)(m + 1) jl - cos C(2n)/(m + 1) ]| ' ■ (^) 

provided 


m > 1 (29a) 

When m = Ij 

) Q = 2 (29b) 

The results of the summation were listed in eguatlons (24) of reference 3 • 
The TTifis.ning of the symbol • E is 

• . S'" 

E = ao® [<£r(l£:i,i - hi^i+x) - «t(ks,i + ks^i+i) I . (30) 

The values of ^i^i+ij ^ 3 , 1 ^ sjod Ics^i+i are obtained from those 

of ki and ks (eq.uatiOES (9a.) and (10a)), respectively, by replacing 
the angle 9 by 2iti/s or 2jt(i + l)/S: 


ki^i = [1 + a cos (2rtni/S)+ (1.6a — 1.8) cos (4nnl/S) 

+ (0.6a — 0.8) cos (6itnl/S) + b sin (2nni/s)- 

+ (3*2b + 3.6ir) sin (4jmi/s) + (l.8b + 2.4jt) sin (6rt7\l/s) ] (3©a) 


» (l/n) [— (2nni/s) — a sin (2jml/S) — (l/2)(1.6a — 1.8) sin (47mi/S) 

— (1/3) (0.6a — 0.8) sin (67tni/S) + b cos (2rtni/s) 

+ ( 1 / 2 ) (3.2b + 3-6rt) cos (4iml/s) + (l/3)(l.8b + 2,hx) cos (6xnl/s)] 

(3Cb) 


In the calculations the E quantities were siumed up numerically. 
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WORK BORE BY THE EXTEEHAL POECES 


It was otserved In. the ea^periments descrihed in reference 9 that 
the stress distrihutlon was not linear in the cutout portion of the cyl- 
inder, although the deviations from linearity were not large as a rule. 

A good appix)ximation to the erparimental curves was obtained by assuming 
a linear force distribution which is not equivalent to a linear stress 
distribution because of the varying amount of effective sheet added to 
the stringer section. Strain-distribution c\irves calculated on the as— 
stanption of a linear force distribution are conpared in figure 4 with 
strains measured in the erperiments . The eapression used for the calcu- 
lation of the force acting upon the ith stringer is 

?i = Per |co8 Ea + (2rti/S) 3 + (S/r) j-^Ecos a + (6/r) ] (31) 

where Per is the c^jc^preesive force acting upon the stringer at the 
edge of the cutout, and 6 is the distance of the neutral axis from the 
horizontal diameter of the cylinder. 

The work done by the external forces acting upon the stringers i« 
equal to the sunajjatlon of the forces times the dleplaceme?».t of the points 
of application of the forces. The displacements of these points are 
equal to the shortening of the distance between the end points of the 
stringers durl^ buckling. Ccnsequently the work is 

L 

W = (1/2) ) Pi r C(6wr/6x)® + (3wt/Sx)^3 dx (32) 

o 

where the summation has to be extended over el 1 the stringers contained 
in one— half the cylinder. Substitutions and integration yield 

L 

W = (l/2)ao^ Pi(rt/L)^ + kg^) P sln^ (2 jcx/l) dx 

= (1/4) (rt2/L)Pcr ao®y (Pi/Pcr) (ki^ + 4^®) (33) 


The summation in the right— hand member of equation (33) was carried out 
numerically , 
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CALCULApCN 0 ? THE BUCSLUfS LOAD 


■ The huolcling ccjndItlorL is 

Ur Ug-trj. + Ustrt; Ufc + Usii = W (3^) 


where the values of the q.uantities mist he taten from eg.uatious (16), 
(l 8 ), (20), (23), (28), and ( 33 )- Eq^uatlon. ( 3 ^) was solved for Per# 
contained in W, hy the use of the following procedure: 


First, a value of n was assumed corresponding to a circumferen- 
tial wave length extending over an integral number of stringer fields. 
With this value M, hi, ks, and were confuted. IText, Per 
assumed. This assumption permitted the calculation of the effective 
width of sheet and ooneeq.uently the moments of inertia of the stringers 
as well as CJeff/Go, The summations were then carried out. Substitution 
of the results in equation ( 3 ^) yields a polynomial of the second degree 
in a and h in the left-hand member, and smother polynomial of the 
second degree in the right-hand member, the latter multiplied by Per* 

Solution for Per gives a fraction wliich can be ^presented symbolically 
as 


Pi(a,b) 
“ P2(a,b) 


( 35 ) 


where pi and P2 stre second— degree polynomials in a and b. The 
values of a eind b, the parameters defining the buckled shape, must 
be chosen so as to make a minimum. It is known from the calculus 

that Per can be minimized by setting 

_ Pi.(a,b) ^ . 

P 2 (a,b) dps/Sa 


The partial differential coefficients of the polynomials pi and ps 
are linear functiens of a and b. Equations (36) r^resent tiiree con- 
nections between Perj b. They were solved by a .rapidly can- 

verging trial-end— error method. First, a and b were calculated from 
the linear equations with the aid of an assumed value of ' Per • Til® val- 
ues of a and b so determined were then substituted into the quadratic 
espression for Per* The procedure was r^eated with the aid of new as— 
sun^tions for Per until the value obtained from the qixadratic eapres— 
Sion was close enough to the assumed value. 
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When the value of obtained in these calculations differed 

materially frcm that aesiimed at the outset, the mcments -of inertia and 
the effective shear modulus had to be calculated again and the entire 
procedure repeated. All the calculations were carried out with different 
values of n. The buchling loads corresponding to these different values 
were conpared, £ind the smallest one was considered as the true buckling 
load. Details of the procedure may be seen from the numerical example 
given in the appendix. 


COMPARISON OF TEEWEY AND EXPERIMENT 


Nvimerical calculations were carried out for the cylinder shown In 
figure 1 for all three of the sizes of the cutout indicated. In each 
case the mlnlmuTn value of the buckling load Per was obtained for 
n = 2 . 666 , that is, when the bulge extended over three stringer fields. 
A typical buckling pattern is shown in figure 2. It .corresponds to a 
cylinder having a 90 ° cutout . 

Some details of the results of the calciilationB are presented in 
table II. The bending moments corresponding to the minimum buckling 
load are plotted against the size of the cutout in figure 5 , which also 
contains the observed bending uShnent at .buckling taken frcm the esperi— 
mental report (reference 9 ) • 


CONCLUSIONS- 


A strain-energy theory has been developed for the calctilation of the 
buckling load in general instability of reinforced circular monocoq.ue 
cylinders which have a symmetric cutout on the compression side and axe 
subjected to pure bending. When -the theory was applied to the teat cyl- 
inders of the earlier experimental report it was found -that at buckling, 
tlie bending moment applied to the cylinders having a 9 ^°, and 135 

cutout was 69-3, 64.8, and 44.5 percent, respectively, of the bending 
moment under which the same type of cylinder buckled when there was no 
cutout. The corresponding -values obtained in the esperimen-bal Investiga- 
tions were 66, 47 ^ and 31 percent when based on the experimental averages, 
and 68 . 2 , 50 « 6 , and 31.3 percent when based on the highest buckling loads 
observed. 


Polytechnic Institute of Brooklyn 
Brooklyn, N. Y., July 8 , 1946. 



NACA TIT No.. 1263 


IT 


APPENDIX 


As a roansrical ezariple of ■buckling load determination, calciilations 
are presented viiicli correspond to cylinders 19 , 20 , and 23 of the test 
series described in reference 9 - The following characteristics of these 
specimens will he needed for calculation: 


Radius, r in. . . 10 

Length of cutout, L.. . in. 19-29 

Distance he'bween adjacent rings, Lx . . . .' in. 6.43 

Wuiaber of rings contained in the length of cutout, m 2 

Ring cross section l/8 hy 3/8 in.. (24s-(I? aluminum, alloj^ 

Number of stringers in full portion of the cylinder, S l6 

Stringer Epaclng along circumference, d in. 3-92? 


Stringer cross section. 
Sheet coTering thickness, t 


in. . . 3/8 X 3/8 (fe4s-^C aluminum alloy) 
in. G.012. . (24S-d? Alclad) 


Cutout angle, 2a deg . 

Young*s modulus, E psi . , IO .5 x 10 

Stringer shear modulus, G . psi .. . 3-9 10' 

Sheet shear modulus at zero ccmqjressive lead, Gq • • psl . . 3-9 ^ 10 


52 


6 


Ccraputations are first given corresponding to an assumed integral 
numher s of stringer fields, say e = 3 , included in the bulge on one 
Bide of the cylinder. The corresponding value of n can be obtained 
from equation (Al) 

n = S/(2s) (Al) 

In this case n = 8/3 = 2.666 . . . Substitution of this value in equa- 
tion (13) gives 


M = 180 , 623.6027 + 3 ^ 8 ^. 5 B 55 a + 93 , 272 . 2882 b. 

+ l,999-4498a2 + 12, l4o . 0843b^ + 2,892.2784ab (A2) 

as can be seen from tab^e I.- Substitution in equation (I 6 ) yields: 

■Or = ( 3 / 16 ) (aoVlOOO) ( 10.5 x 10® x 80.35 x 10~®) (2 + l)M (A3) 

where the moment of inertia I^ of the ring cross sectiem augmented by 
an effective width of sheet (taken equal to the width of the ring) is 


(a4) 


^r = ,( 1 / 12 ) ( 0.125 + 0 . 012 )® (0.375) = 80.35 X 10 “® in.^ 



18 


MCA TK NO. 1263 


With, the aid of e^iiatlon (A2), the ring strain energy hecoBies: 

Ur = ao^ ( 85 , 718. 0013 + 4 , 190 . 2299 a + 44 , 263 . 9 ^ 00 b 

+ 9^.8729a2 + 5 , 761 . 2834h2 + 1,372. 5799a'b) . (A 5 ) 


The functions ki, ks, k 4 can he calculated from, the assumed value 

of n. . Nor this purpose it is convenient to arraiage the trigonometric 
functions needed and their coefficients in tabular form, as shown in 
table III. 


In this table the first three rows contain the trigonometric func- 
tions of the angle .cp. 'As only values of k^, kg, ki which correspond 
to integreil numbers of stringer fields will be used in the summations to 
be evaluated later, the angle <p was replaced by its equivalent (2ni/s) , 
where 2«/s is the angle subtended by one stringer field. 

The third, fourth, and fifth rows contain the coefficients of the 
trigonometric functions appeegrlng in each column above them. These coef- 
ficients are different for kx, ka, and k 4 , and can be obtained fsrbm 
equations (9a), (lOa), and (22a), respectively. 

The value of kx for i = 0 Is obtained by multiplying the ezpres— 
eions appearing in the same column in the first and fourth rows, aaad by 
adding the resulting eight products. Products of the elements of the 
second and fourth rows will lead to the value of kx for i = 1 , and 
similarly for all others. The results are tabulated in table IV. 

The next step is the aseuuqitlon of the critical load . for the 

purpose of obtaining the effective width of sheet 2 w to be added to 
the stringers, the moments of inertia of the stringers, the effective 
shear modulxxs &eff , shift of the neutral axis from the horizoir- 

tal axis of the cylinder 6 . From an assumed value of P^x. = 3370-5 
8 was found to be 2.4 in.; the other quantities are listed in table V. 

In table V column (2) Is obtained from equation (31), column. (5) by di- 
viding column ( 2 ) by E = IO .5 X 10® psl. Columns (4), ( 5 ), and ( 6 ) can 
be most conveniently obtained by the use of a previously drawn curve of 
the strain c against the area Aeff of stringer and effective width 
combination. A ciirve of this type was used in the present calculations 
6 ind was constructed with the aid of the foUow.ing formula for effective 
width: 


2w =.(l/e)(d/r) -|p.3t -h 1.535 


(t/d) (er — 0 . 3 t)r 


1 / 2 ' 
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Coliunns (8) were calculated according to the principles stated in. 

reference 1 . In. column ( 9 ) the te Iu^ of the budcling load of the sheet 
Spy was calculated to he 3*3 '• 10 • Coluam. (lO) was obtained from 

eijuation (25) . 


The strain energy stored in the stringers because of bending in the 
radial direction is from equation (I8) : 


U. 


r 


str^, = ^o^ /(I9.29) 


10.5 X 10 


I 

i=o 


I, 


StTr 


(A6) 


since the modulus E is constant for all stringers. The summation in- 
dicated is the sum of the product of the values in column (7) in table V 
for any value of i, and the- squares of the quantity ki in table IV 
for the same value of 1. ^^etrt obtained by substituting coluam (8) 

for column (7) in table V and kg for ki in table IV. These opera- 
tions yield: 

^str = Ustrr + ^str^ = ao^ (102, 882. 5969 - 1,453.3391a 

+ 53,374.568Cfb + 3,945. 3l60a^ + 7,242. 3553b^ - 409.40^ab) . (A7) 


The torsional strain energy in the stringera contains the sl frmn atlcn 

as can be seen from equation ( 23 ) . This is the sum of the squares 

of the values of k^ given in table XV. The complete result for the 
torsional strain energy wae found to be: 

o't = 80^(300,981.0280 + 25,092.1760a + 1^, 584.300b 

+ I,656.71l6a2 +. 22 ,a 30 , 3894 b^ +. 6 , 426 : 6 l 0 ab) (AS) 


The expression given in equation ( 28 ). for the strain energy due to 
shear in the sheet covering required the evaluation of the quantity -Ei. 
This quantity, .given in .equation (30), can be easily obtained if it is 
noticed that the terms hi,l ~ and ks^i + ks^i+i are respec- 

tively the difference the simi of tenns appearing In adjacent rows of 
the kj. and kg columns of table IV. ¥ith the values of Geff/Oo 
taken from table Y and with = 0.03927 and c&fc = ” 0.5 from equation 
(27) the shear strain energy was found to be: 
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■f^sh “ ao^(l6,6k3 .k07& + 6,t(l6.0112a + 5,515 •529fj'b + l,l80.0227a^ 

+ 1,232 .0806-b^ + l88.3523a‘b) (A9) 

The work done hy the extonial forcea is given in equation (33) • Sach 
terra of the aumiaation contaluad in that equation la the aum of the aquarea 
of the vaLuoa of k and k (ta'Dle 17) multiplied Tby which 

can he ob-oained from table Y. The reault of theae calculationa was: 

W = Per 80^(17*522403 - 2.59it-795a + 9*796509t 

+ 1.422668a^ + 1.426l08b^ - 0.429835ab) (AlO) 


The value a of the strain ener^ea and external work token from equa- 
tions (A5), (A7), ,(a8), (A9), and (AlO) were substituted in equation (34), 
with the reault: 

ao"(506,227.034 + 3^, 244.578a + 268,738.347t 
+ 7,730*9227h^ + 37,066.1019b® + 7,578.13570b) 

= Oq® Poy( 17. 522403 - 2.594795a + 9.796509b 

+ 1.422668a® + 1.426108b® - 0.4298350b) (All) 

According to equation (36), equation (All) was solved for Pcr> '4bo 

numberator and denominator of the resulting expression were differenti- 
.ated with respect to a and b. The restilt is given in equations (A12) : 




3!).. 244.578 + 15,^161. 8454a + 7.578.13' 

-2.594795 + 2.845336a - 0.4^6351 


268,738.347 t 7.578.13570 t 74.132.2038b 
9.796509 - 0.429835a + 2.852216b 


■^cr 


> 


(A12) 


cr 


Theae equations were reduced to two linear equations in a and b by 
assuming a value of, Per ® 3770 and clearing fractions. These equations 
were solved simultaneously for a and b, with the following result: 


a = -3.0541 


(AI3) 


b » -3.2142 
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Substitution of these values In equation (All) led to a new value of 
I’cr = 3793 Ih. This restilt was considered to he sufficiently close to 
the original assumption to make it unnecessary to repeat the calculations 
for this value of n. 

Repetition of the entire procedure for values of n corresponding 
to s = 2 and s = 4 gave Per “ 11544 lb and = 3902 lb, respec- 
tively. The value Per = 3793 lb was the lowest of the three and was 
therefore considered the true buckling load. 
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TABEE I.- VATOBS OF M 


8 

n 

Constsuxt 

a 

1 

a 

a 


ab 

i 

k 

624952,2978 

31194.1304 

322964.0316 

6969.3405 

42066.8361 

10160.6399 

B 

2.666... 

IB0623.6027 

8829.5855 

93272.2882 

1999.4498 

12140.0843 


B 

2 

73562.5088 

3487,1679 

37945.9796 

805.8166 

4933.8711 


5 

1.6 

35969.3603 

1634.3092 

18527.8397 

388.6488 

2405.8392 

546 . 4166 .. 

6 

1 . 333 ».. 

19654.0645 

843.6900 

10105.7565 

208.7710 

1310,0700 

286,81^; 
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TABLE n 


•p- 


1 ^ 

Size of 
cutout 
(deg) 

2 n 

n 

a 


^cr 

*^or 

calculated 

Her 

avenage 

expertraental 

Peroent 
difference 
of eaeperi- 
mental 
average 

M^st 

eerperi— 

mental 

moment 

1 

Percent 
difference 
of highest 
experimental 
moment 

^5 

2 

4 

o- 3 on 

iPi 

13,564 







3 

2 > 666 ... 

-3.1734 


3,722 

258,109 

192,800 

33.9 

197,600 

30.6 


4 

2 

-2.7612 

- 3 -S 710 

3,875 






90 

2 

4 

.3417 

-3.6809 

U ,544 







3 1 

2,666-,., 

-3.0541 

- 3.2142 

3,793 

187,216 

135, SQO 

37.9 

146,000 

28.2 


4 

2 

-2.724s 

- 3.2831 

3,902 






135 - 

2 

4 

*3921 

-3.6881 

14,549 







3 

2.666*,. 

-2.9556 

-3.2553 

3,847 

120,471 

89,500 

43-5 

92,200 

39.3 


4 

2 

-2.6519 

-3.3001 

3,956 
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TABLE IV 








! 0 
1 

- 1.6 + 3 . 2 a 

3.0630528 + 1.2B 

120.637166 + 34. 13333.. .b 

i 

j 

12.4945177 “ 0.9a + 3.6373O7IB 

-2.10318072 - 0.5845672a - O.3375B 

- 82.164026 - 9.699486a - 21.6b 

i2 

1 

t 

-6.6945177 - 0.7a - 1,905256113 

-1.19549133 - 0.06495191a - 0.2605b 

21.845445 + 5.080683a + 4 . 533 ...b 


TABLE V 


jd) 

i 

( 2 ) 

Pi 

( 3 ) 

€ A 

( 4 ) 

e 

( 5 ) 

Aeff 

( 6 ) 

2 w 

( 7 ) 

Pstrj. 

( 8 ) 

^strt 

( 9 ) 

( 10 ) 

0 

3370.5 

3.21 X 10 "* 

21.4 X 10 *^ 

0.1505 

1.567 

2025 X 10® 

■ 3 , 500 - X 10 ® 

6.06 

0,495 

1 

2215.815 

2.1103 

12.8 

■ •1^. 

2 . OCX) 

2340 

10,000 

3 . 88 - 

.562 

2 

853.75 

0.8131 

4.33 

.I&78 

3.927 

2392 • 

62,200 

1,312 

.745 
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SECTION A-A COMMITTEE FD» AERONAUTICS 

FIGURE I- MONOCOQUE CYLINDER. 
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Pig. 2 


n =? 2.666... 

a = “3 .0541 
b = -3 .2142 
S = 3 
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FIGURE 2.- TYPICAL RADIAL DEFLECTION PATTERN. 

DEFLECTION EXAGGERATED 
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Fig. 4 



FIGURE 4 -COMPARISON OF EXPERIMENTAL AND 

ASSUMED STRAIN DISTRIBUTION/ 
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